Substituted 3-phenylpropenoates and related analogs: electron ionization mass spectral fragmentation and density functional theory calculations.

Analysis of ethyl 3-(2-chlorophenyl)propenoate by electron ionization mass spectrometry showed the distinct loss of an ortho chlorine. To characterize the structural requisites for the observed mass fragmentation, a series of 30 halogen-substituted 3-phenylpropenoate-related structures were examined. All ester-containing alkene derivatives exhibited loss of the distinctive chlorine from the 2-position of the phenyl ring. Analogous derivatives with the halogen (chlorine or bromine) in the para position did not evidence selective halogen loss. Results demonstrated that substituted 3-phenylpropenoates and their analogs fragment via the formation of a previously reported benzopyrylium intermediate. To understand the correlation between the intramolecular radical substitution and the abundance and selectivity of the chlorine (or other halogen) displacement, density functional theory calculations were performed to determine the charge on the principal cation involved in the chlorine loss (in the ortho, meta, and para positions), the charge for the neutral radical (noncation), the excess alpha-electron density on the relevant atom and the energy to form the cation from the neutral atom (ionization energy). Results showed that the selectivity and extent of halogen displacement correlated highly to the electrophilicity of the radical cation as well as the neutral radical. These data further support the proposed fragmentation mechanism involving intramolecular radical elimination.

A QSAR for the mutagenic potencies of twelve 2-amino-trimethylimidazopyridine isomers: structural, quantum chemical, and hydropathic factors.

An isomeric series of heterocyclic amines related to one found in heated muscle meats was investigated for properties that predict their measured mutagenic potency. Eleven of the 12 possible 2-amino-trimethylimidazopyridine (TMIP) isomers were tested for mutagenic potency in the Ames/Salmonella test with bacterial strain TA98, and resulted in a 600-fold range in potency. Structural, quantum chemical, and hydropathic data were calculated on the parent molecules and the corresponding nitrenium ions of all of the tested isomers to establish models for predicting the potency of the unknown isomer. The principal determinants of higher mutagenic potency in these amines are: (1) a small dipole moment, (2) the combination of b-face ring fusion and N3-methyl group, (3) a lower calculated energy of the pi electron system, (4) a smaller energy gap between the amine HOMO and LUMO orbitals (Pearson "softness"), and (5) a more stable nitrenium ion. Based on predicted potency from the average of six regression models, the isomer not yet synthesized and tested is expected to have a mutagenic potency of 0.77 revertants/microg in tester strain TA98, which is near the low end of the potency range of the isomers.

Experimental and simulation studies of heat flow and heterocyclic amine mutagen/carcinogen formation in pan-fried meat patties.

Heterocylic amine (HA) compounds formed in the cooking of certain foods have been shown to be bacterial mutagens and animal carcinogens, and may be a risk factor for human cancer. To help explain the variation observed in HA formation under different cooking conditions, we have performed heat-flow simulations and experiments on the pan-frying of beef patties. The simulations involve modeling the heat flow within a meat patty using empirically derived thermal transport coefficients for the meat. The predicted temperature profiles are used to integrate the Arrhenius rate equation to estimate the concentration of HAs formed in the meat. We find that our simulations accurately model experimentally determined temperature profiles, cooking times, HA spatial distributions and total HA formation in patties that are flipped once during the pan-frying process. For patties flipped every 60 s, the simulations qualitatively agree with experiment in predicting reduced cooking times and HA formation relative to the singly-flipped patties. However, the simulations overestimate the effect of rapid flipping on cooking times and underestimate the effect of flipping on total HAs formed. These results suggest that the dramatic reductions in HA formation due to rapid flipping may be due to factors other than the heating process or that there is a critical feature of the flipping process that is not captured in our model.

Immunoreactivity of organic mimeotopes of the E2 component of pyruvate dehydrogenase: connecting xenobiotics with primary biliary cirrhosis.

In primary biliary cirrhosis (PBC), the major autoepitope recognized by both T and B cells is the inner lipoyl domain of the E2 component of pyruvate dehydrogenase. To address the hypothesis that PBC is induced by xenobiotic exposure, we took advantage of ab initio quantum chemistry and synthesized the inner lipoyl domain of E2 component of pyruvate dehydrogenase, replacing the lipoic acid moiety with synthetic structures designed to mimic a xenobiotically modified lipoyl hapten, and we quantitated the reactivity of these structures with sera from PBC patients. Interestingly, antimitochondrial Abs from all seropositive patients with PBC, but no controls, reacted against 3 of the 18 organic modified autoepitopes significantly better than to the native domain. By structural analysis, the features that correlated with autoantibody binding included synthetic domain peptides with a halide or methyl halide in the meta or para position containing no strong hydrogen bond accepting groups on the phenyl ring of the lysine substituents, and synthetic domain peptides with a relatively low rotation barrier about the linkage bond. Many chemicals including pharmaceuticals and household detergents have the potential to form such halogenated derivatives as metabolites. These data reflect the first time that an organic compound has been shown to serve as a mimeotope for an autoantigen and further provide evidence for a potential mechanism by which environmental organic compounds may cause PBC.

Extended quantitative structure-activity relationships for 80 aromatic and heterocyclic amines: structural, electronic, and hydropathic factors affecting mutagenic potency.

The mutagenic/carcinogenic heterocyclic amines formed during the cooking of protein foods have been determined to be probable or possible human carcinogens. As part of a comprehensive study of the food mutagens, our laboratory has produced a series of quantitative structure-activity relationships (QSARs) of aromatic and heterocyclic amines, to attempt to elucidate the mechanisms of mutagenesis/carcinogenesis. Amines are genotoxically active only after activation by a series of reactions converting the parent compound to an electrophilic derivative, which is postulated to be a nitrenium ion that covalently binds to and damages DNA. An important agent in this conversion is cytochrome P450. In this report we develop a QSAR for 80 amines of diverse structure and a range of 10 orders of magnitude in mutagenic potency. New structural factors and quantum chemical ab initio and Hückel calculations are included. The results are interpreted to show that a main determinant of mutagenic potency is the extent of the aromatic pi-electron system. Small contributions are made by both the dipole moment and the calculated stability of the nitrenium ion. Multiple linear regression models account for nearly two-thirds of the variance in potency, leaving room for additional unknown factors. The role of cytochrome P450 1A in amine toxification is supported, and further theoretical and experimental research on its reaction mechanisms and modeling of its active site are proposed.

Quantitative structure-activity relationship of flavonoids for inhibition of heterocyclic amine mutagenicity.

The mutagenic/carcinogenic heterocyclic amines formed during the cooking of protein foods have been determined to be a potential risk to human health. Therefore, mitigation measures are beginning to be studied. A recent finding is that the induction of mutation in Salmonella by these amines can be inhibited by the addition of flavonoids to the assay. This study combines data on the inhibitory process with structural, ab initio quantum chemical, hydropathic, and antioxidant factors to develop a quantitative structure-activity relationship (QSAR) database and statistical analysis. For 39 diverse flavonoids the inhibitory potency varied approximately 100-fold. Three predictive variables, in order of decreasing contribution to variance, are: (1) a large dipole moment; (2) after geometric minimization of energy, a small departure from planarity (i.e., small dihedral angle between the benzopyran nucleus and the attached phenyl ring), and a low rotational energy barrier to achieving planarity; and (3) fewer hydroxyl groups on the phenyl ring. However, these variables account for less than half of the variance in inhibitory potency of the flavonoids. Frontier orbital energies and antioxidant or radical scavenging properties showed little or no relationship to potency. We conclude that interference by the flavonoids with cytochrome P450 activation of the promutagens is the probable mechanism for inhibition of mutagenesis, and suggest avenues for further research. Environ. Mol. Mutagen. 35:279-299, 2000 Published 2000 Wiley-Liss, Inc.

Heterocyclic amine formation and the impact of structure on their mutagenicity.

The occurrence and formation of heterocyclic amines in foods is discussed in light of the consistent finding of a new class of imidazopyridines. In addition, a quantitative structure-activity relationship will be presented correlating the potency of these imidazopyridines to predicted chemical properties. Although no strong linear correlation is found between the potency and the chemical properties, a low dipole moment is found to be a qualitative predictor of high mutagenic potency.

Chemical factors in the action of phosphoramidic mustard alkylating anticancer drugs: roles for computational chemistry.

The nitrogen mustard based DNA alkylating agents were the first effective anticancer agents and remain important drugs against many forms of cancer. More than fifty years of research on the nitrogen mustards has yielded a broad range of therapeutically useful compounds and a detailed knowledge of the biochemical mechanism of these drugs. Nevertheless, there is much ongoing research on the phosphosphoramidic and other nitrogen mustards to increase their potency and reduce their toxic and mutagenic side effects. To understand the existing nitrogen mustards, and to design the next generation of these drugs, more knowledge is needed about the effects of chemical modifications on their activation and selectivity. Because of the existing knowledge of these drugs, atomic-level chemical modeling can play an important role in the understanding of the phosphoramidic mustard compounds; however, it has not proved straight forward to directly relate the activity of these mustards with simple chemical properties such as bond lengths or atomic charges. Instead, quantum chemical simulations will be required to simulate the activation and alkylation reactions of these compounds, which will require the newest generation of quantum chemical and solvent modeling methods. Additionally, molecular dynamics simulations of the adducted DNA can provide data on the factors favoring crosslinking and its structural consequences. This review summarizes the extensive literature on the metabolism, activation, and action of the phosphoramidic mustards, with an emphasis on the roles that chemical modeling has and will play in the development of this important class of drugs.

Interaction and solvation energies of nonpolar DNA base analogues and their role in polymerase insertion fidelity.

Although DNA polymerase fidelity has been mainly ascribed to Watson-Crick hydrogen bonds, two nonpolar isosteres for thymine (T) and adenine (A)--difluorotoluene (F) and benzimidazole (Z) --effectively mimic their natural counterparts in polymerization experiments with pol I (KF exo-) [JC Morales and ET Kool. Nature Struct Biol, 5, 950-954, 1998]. By ab initio quantum chemical gas phase methods (HF/6-31G* and MP2/6-31G**) and a solvent phase method (CPCM-HF/6-31G**), we find that the A-F interaction energy is 1/3 the A-T interaction energy in the gas phase and unstable in the solvent phase. The F-Z and T-Z interactions are very weak and T-Z is quite unstable in the solvent. Electrostatic solvation energy calculations on F, Z and toluene yield that Z is two times, and F and toluene are five times, less hydrophilic than the natural bases. Of the new "base-pairs" (F-Z, T-Z, and F-A), only F-A formed an A-T-like arrangement in unconstrained optimizations. F-Z and T-Z do not freely form planar arrangements, and constrained optimizations show that large amounts of energy are required to make these pairs fit the exact A-T geometry, suggesting that the polymerase does not require all bases to conform to the exact A-T geometry. We discuss a model for polymerase/nucleotide binding energies and investigate the forces and conformational range involved in the polymerase geometrical selection.

Chemical and biological factors affecting mutagen potency.

This review surveys the chemical and biological factors that are correlated with the mutagenic activity of the aromatic and heterocyclic amines. Particular attention is given to the predicted quantum chemical properties of the parent amines and their metabolites. A number of chemical properties have been found to correlate well with measured mutagenic potency, including log P, the energies of the frontier orbitals of the parent amines, and the thermodynamic stability of the nitrenium ion, possibly the ultimate DNA-binding species. These correlations are intriguing clues to the mutagenic activity of the aromatic amines; however, many factors still await final explanation, including the exact mechanisms of the metabolic enzymes, the identity(s) of the ultimate DNA-binding species, the reaction mechanism in the DNA-adduction, the role of sequence context in the covalent and non-covalent binding of the adducts, and the role of DNA repair.

Elements in abasic site recognition by the major human and Escherichia coli apurinic/apyrimidinic endonucleases.

Sites of base loss in DNA arise spontaneously, are induced by damaging agents or are generated by DNA glycosylases. Repair of these potentially mutagenic or lethal lesions is carried out by apurinic/apyrimidinic (AP) endonucleases. To test current models of AP site recognition, we examined the effects of site-specific DNA structural modifications and an F266A mutation on incision and protein-DNA complex formation by the major human AP endonuclease, Ape. Changing the ring component of the abasic site from a neutral tetrahydrofuran (F) to a positively charged pyrrolidine had only a 4-fold effect on the binding capacity of Ape. A non-polar 4-methylindole base analog opposite F had a <2-fold effect on the incision activity of Ape and the human protein was unable to incise or specifically bind 'bulged' DNA substrates. Mutant Ape F266A protein complexed with F-containing DNA with only a 6-fold reduced affinity relative to wild-type protein. Similar studies are described using Escherichia coli AP endonucleases, exonuclease III and endonuclease IV. The results, in combination with previous findings, indicate that the ring structure of an AP site, the base opposite an AP site, the conformation of AP-DNA prior to protein binding and the F266 residue of Ape are not critical elements in targeted recognition by AP endonucleases.

Quantitative structure-activity (QSAR) relationships of mutagenic aromatic and heterocyclic amines.

We extended our previous studies of mutagenic/carcinogenic heterocyclic aromatic amines formed during the cooking of foods to 66 aromatic and 99 heterocyclic amines for which mutagenic potency data are available. The amines require activation by enzymes to form metabolites reactive with DNA and exhibit an enormous range of potency as frameshift mutagens in the Ames/Salmonella assay. To ascertain factors that might influence potency, structural features and quantum mechanical parameters calculated by the Hückel method (and, for a subset of 20 amines, by semi-empirical AM1, and ab initio methods) were analyzed by multiple linear regression. The major findings were: (1) earlier findings on cooked food mutagens and their synthetic congeners can be extended to other amines; (2) mutagenic potency is directly related to the number of fused aromatic rings (size of the aromatic system), the number of ring nitrogen atoms (participation of lone electron pairs in the pi-cloud), and presence of a methyl substituent on a ring nitrogen; (3) potency is inversely related to the energy of the lowest unoccupied molecular orbital (LUMO) of the parent amine. Ford and Griffin (1992) and Sabbioni and Wild (1992) showed that the LUMO energy of the derived nitrenium ion is closely related to its stability (calculated with reference to aniline). Increased stability has been hypothesized to enhance the probability of adduct formation with DNA by avoiding detoxifying side reactions and increasing the lifetime of the ion. In the large heterogeneous series of amines in our present study the Hückel method energy of the highest occupied molecular orbital (HOMO), rather than the LUMO energy, of the nitrenium ion was marginally related to the potency of the parent amine. However, in the selected subset of 20 amines with ab initio calculation, the LUMO energy of the ion confirmed the previous reports. The contribution of quantum chemical factors to mechanistic insight on the mutagenicity and carcinogenicity of aromatic and heterocyclic amines is still under development.

Hydration effects on the duplex stability of phosphoramidate DNA-RNA oligomers.

Recent studies on uniformly modified oligonucleotides containing 3'-NHP(O)(O-)O-5'internucleoside linkages (3'amidate) and alternatively modified oligonucleotides containing 3'-O(O-)(O)PNH-5'internucleoside linkages (5'amidate) have shown that 3'amidate duplexes, formed with DNA or RNA complementary strands, are more stable in water than those of the corresponding phosphodiesters. In contrast, 5'amidates do not form duplexes at all. There is no steric reason that the 5'amidate duplex should not form. We demonstrate that these differences arise from differential solvation of the sugar-phosphate backbones. By molecular dynamics calculations on models of 10mer single-stranded DNA and double-stranded DNA-RNA molecules, both with and without the phosphoramidate backbone modifications, we show that the single-stranded 3'amidate and 5'amidate backbones are equally well solvated, but the 5'amidate backbone is not adequately solvated in an A-form duplex. These results are supported by quantum chemical free energy of solvation calculations which show that the 3'amidate backbone is favored relative to the 5'amidate backbone.

Deprotonation and hydride shifts in nitrenium and iminium forms of aminoimidazole-azaarene mutagens.

The mutagenicity of many 2-aminoimidazole-azaarenes (AIA) is thought to be mediated by the nitrenium form of the exocyclic amine. This hypothesis is supported by the numerous correlations found between calculated and experimentally-measured chemical properties for the nitreniums and the mutagenic potencies of the nitreniums and their parent amines. One factor favoring high mutagenic potency is the presence of a methyl substituent in the 1- or 3-imidazole position. In this paper, we investigate both the deprotonation of the imidazole ring nitrogens in non-N-methylated AIA mutagens and the plausibility of a chemical pathway involving a 1-4 hydride shift to form an iminium ion, thereby stabilizing the cationic N-methyl substituted AIA mutagens. It has been widely noted that factors that stabilize the nitrenium moiety lead to significantly higher mutagenic potency; hence, the transformation of the nitrenium to a more stable species might be expected to increase the potency, provided that it does not eliminate the electrophilic reactivity of the compound. Using ab initio quantum chemistry and polarizable continuum solvation models, we find that the imidazole ring nitrogens of the nitrenium ions are extremely acidic. This suggests that upon formation of the exocyclic nitrenium these sites will deprotonate to form a neutral imine. We have also studied the 1-4 hydride shift from an imidazole ring methyl to the exocyclic nitrenium to form an iminium. We predict that for AIA mutagens with just two fused rings the resulting iminium species are more stable in the gas phase than the corresponding nitreniums. For mutagens with larger conjugated systems, the nitrenium is stabilized by resonance and is more stable than the corresponding iminium. In the aqueous phase, however, the iminium form is predicted to be more stable than the nitreniums for all polycyclic compounds studied. Although equilibrium calculations favor the iminium form, these have been experimentally shown to be short-lived and their actual concentration will depend on the complex kinetics of AIA mutagen metabolism. The quantum chemical results also show a strong correlation between the relative iminium-nitrenium energy difference and the charge on the exocyclic nitrogen.

Comparison of crystal structure and theory for 2-amino-1-methyl-6-phenylimidazo[4,5-b]pyridine.

The crystal structure of the food mutagen 2-amino-1-methyl-6-phenylimidazo[4,5-b] pyridine (PhIP) has been determined by single-crystal X-ray crystallography. Crystals grown by evaporation of an aqueous solution form in the monoclinic space group P2(1)/n with two molecules of PhIP per asymmetric unit, along with six water molecules. The phenyl groups of these two PhIP molecules have torsion angles of different magnitude with respect to the plane of the imidazopyridine moiety. To maintain centrosymmetry, the crystal also contains an oppositely torsioned symmetry equivalent of each. The amino groups of both PhIP molecules take part in an extensive hydrogen bond network with the water of crystallization, forming long channels through the crystals parallel to the crystallographic b axis. The diffraction results are compared to theoretical calculations of the optimized geometry for a single PhIP molecule in vacuo as well as with water hydrogen-bonded to the exocyclic amine. In general, the agreement between the X-ray crystal structure of PhIP and its theory-derived counterpart in vacuo is within the combined experimental-theoretical uncertainty. The C-N bond to the exocyclic amine and the neighboring C=N imidazole bond are exceptions. This is attributed to the combined neglect of the crystal environment, waters of hydration, and the lack of coplanarity between the imidazole ring and the amine group in the calculations. To address the effect of waters of hydration, additional calculations were performed to optimize the geometry of a PhIP molecule with two water molecules hydrogen-bonded to the exocyclic amine. The resulting C-N exocyclic amine and C=N imidazole bond lengths were closer to those obtained by X-ray diffraction. The accord between theory and experiment demonstrates the utility of applying theory to (1) accurately predict structures of PhIP metabolites and intermediates that are too labile for study by conventional structural techniques such as X-ray crystallography and (2) assist in studying the mechanisms by which PhIP and its metabolites interact with proteins and DNA.

Structural and quantum chemical factors affecting mutagenic potency of aminoimidazo-azaarenes.

A set of 16 mutagenic aminoimidazo-azaarenes, including four that have been isolated from cooked foods and identified as bacterial mutagens and rodent carcinogens, was selected from a larger series previously published [Hatch et al. (1991): Environ Mol Mutagen 17:4-19] for an in-depth structure-activity study using computational methods. Structural features believed to affect mutagenic potency were tabulated. Molecular orbital energies and other electronic properties of these compounds were calculated using Huckel, semiempirical AM1, and ab initio quantum mechanical methods. Factor interrelationships were studied by multiple linear regression and canonical correlation analyses. Our goal was an improved understanding of the chemical basis of mutagenicity for this class of heterocyclic amines. The major findings were as follows: 1) mutagenic potency is related to the size of the aromatic ring system; 2) potency is enhanced by the presence and location of an N-methyl group; 3) potency is enhanced by addition of ring nitrogen atoms in pyridine, quinoline, and quinoxaline configurations; 4) potency is inversely related to the energy of the LUMO (lowest unoccupied molecular orbital) of the parent amines; 5) potency is directly, though weakly, related to the LUMO energy of the derived nitrenium ions; and 6) the calculated thermodynamic stability of the nitrenium ions (relative to the parent amine) is directly correlated with nitrenium LUMO energy and with the negative charge on the exocyclic nitrogen atom. Although this study raises several intriguing issues relating mutagenicity to chemical properties, further study will be required to determine the plausibility of the nitrenium ion as the ultimate mutagen for binding to DNA.

A structural basis for a phosphoramide mustard-induced DNA interstrand cross-link at 5'-d(GAC).

Phosphoramide mustard-induced DNA interstrand cross-links were studied both in vitro and by computer simulation. The local determinants for the formation of phosphoramide mustard-induced DNA interstrand cross-links were defined by using different pairs of synthetic oligonucleotide duplexes, each of which contained a single potentially cross-linkable site. Phosphoramide mustard was found to cross-link dG to dG at a 5'-d(GAC)-3'. The structural basis for the formation of this 1,3 cross-link was studied by molecular dynamics and quantum chemistry. Molecular dynamics indicated that the geometrical proximity of the binding sites also favored a 1,3 dG-to-dG linkage over a 1,2 dG-to-dG linkage in a 5'-d(GCC)-3' sequence. While the enthalpies of 1,2 and 1,3 mustard cross-linked DNA were found to be very close, a 1,3 structure was more flexible and may therefore be in a considerably higher entropic state.

Comparison of the protonation of isophosphoramide mustard and phosphoramide mustard.

The alkylating agent isophosphoramide mustard (IPM) spontaneously forms a relatively stable aziridine derivative which can be directly observed using NMR spectroscopy. The protonations of IMP and its aziridine were probed using 1H, 31P, 15N, and 17O NMR spectroscopy. The positions of the 31P, 15N, and 17O resonances of IPM between pH 2 and 10 each exhibit a single monobasic titration curve with the same pKa of 4.31 +/- 0.02. On the basis of a comparison with other compounds and our earlier work with phosphoramide mustard, the NMR results for IPM indicate that protonation occurs at nitrogen and not oxygen. Over this same pH range, each of the 1H, 31P, and 15N resonances of IPM-aziridine also show a single monobasic titration with a pKa of 5.30 +/- 0.09. The magnitude of the change in chemical shifts suggests that the protonation of the IPM-aziridine occurs at the ring nitrogen. Theoretical gas-phase calculations of PM, IPM, and IPM-aziridine suggest O-protonation to be more likely; however, aqueous phase calculations predict the N-protonated forms to be most stable. Furthermore, for PM and IPM-aziridine, which contain nonequivalent nitrogens, the theoretical calculations and experimental data both agree as to which nitrogen undergoes protonation. These results suggest that the IMP-aziridine remains unprotonated under physiological conditions and may, in part, explain the lower alkylating activity of IPM as compared to PM.

Molecular pharmacology of hepsulfam, NSC 3296801: identification of alkylated nucleosides, alkylation site, and site of DNA cross-linking.

We have determined that hepsulfam, in common with its structural homologue busulfan, alkylates both free guanosine and GMP in DNA at the 7 nitrogen. Mass spectral analysis of the products of the reaction of hepsulfam with guanosine has identified the mono- and bis-alkylated guanosine adducts. UV spectrophotometry and mass spectrometry were used to confirm that alkylation occurred at the 7 nitrogen by following the formation of the formamidopyrimidyl form of the hepsulfam-guanosine adduct at high pH. We have also isolated and identified 1-guanyl,7-hydroxyheptane, 1-guanyl,7-sulfamylheptane, and 1,7-bis(guanyl)heptane from in vitro reaction mixtures of hepsulfam and calf thymus DNA. We have isolated bis-(7-formamidopyrimidyldeoxyguanosinyl)-heptane from an enzymatic digest of DNA treated with hepsulfam. Finally, we have found that hepsulfam forms interstrand cross-links at 5'-GXC-3' sites in model oligonucleotides.

Protonation of phosphoramide mustard and other phosphoramides.

The chemistry of the bifunctional alkylating agent phosphoramide mustard and model phosphoramides was probed by multinuclear NMR spectroscopy as a function of pH. Between pH 1 and 11, both the 31P and 15N resonances for phosphoramide mustard displayed a single monobasic titration curve with a pKa of 4.9. The protonation below pH 4.9 correlates with the loss in reactivity of the mustard. The 17O NMR spectrum of 17O-enriched phosphoramide mustard shows little change with pH. The data on the mustard was compared to 15N and 31P NMR data on 15N-enriched phosphoramidic acid, phosphorodiamidic acid, and phosphoric triamide. Contrary to the conclusions of previous studies, our combined 31P, 15N, and 17O NMR results are more consistent with N-protonation of phosphoramide mustard rather than an O-protonation. Theoretical calculations on the phosphoramidic acid, phosphorodiamidic acid, and phosphoric triamide show O-protonation to be more stable in the gas phase. For the latter two compounds, the calculations suggest that N-protonation may be the most stable protonated form in the aqueous phase. These findings influence our understanding of the structure-activity relationships of phosphoramide mustards.